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We consider the detailed structure of low energy excitations in the periodic spin-1/2 XXZ Heisen- 
berg chain. By performing a perturbative calculation of the non-linear corrections to the Gaussian 
model, we determine the exact coefficients of asymptotic expansions in inverse powers of the system 
length N for a large number of low-lying excited energy levels. This allows us to calculate eigenen- 
ergies of the lattice model up to order 0(N~ 4 ), without having to solve the Bethe Ansatz equations. 
At the same time, it is possible to express the exact eigenstates of the lattice model in terms of 
bosonic modes. 

PACS numbers: 75.10.Pq, 05.30-d, 02.30 Ik 

I. INTRODUCTION 

Quantum models defined on discrete lattices are very common in solid state theory. Two routes to study their 
properties are conceivable: Either the attempt to solve the model on the lattice, or the formulation of an effective field 
theory in the continuum. Lattice models arc much more tractable than continuum theories in numerical simulations. 
On the other hand, within the field-theoretical picture, it is often possible to describe the spectrum in terms of non- 
or weakly interacting quasiparticles, which makes this approach very attractive from an analytical point of view. 
Thus, it is most desirable to express the lattice eigenstates in terms of the conceptually much simpler field theoretical 
eigenstates. 

This goal is generally not achievable, since neither the lattice model nor the full field theory can be solved without 
approximations. Integrable one-dimensional models, however, are the most promising candidates where such a de- 
scription can be realized quantitatively. Indeed, the effective bosonic theory of the spin-1/2 XXZ Heisenberg chain 
offers the opportunity to obtain the lattice eigenenergies as an asymptotic expansion in the inverse system length, 
and at the same time to express the exact lattice eigenstates as linear combinations of bosons. So far, finite-size cor- 
rections to the bosonic spectrum have been determined quantitatively for particle excited states. 1,2 We now calculate 
the coefficients of the asymptotic corrections quantitatively for current and particle-hole excited states as well. 

In particular, we obtain the leading terms in an expansion of the lattice energies in the inverse system length to 
order O (./V -4 ) , which translates into a relative deviation of a fraction of a percent or less already for moderate chain 
lengths N ~ 20 — 100. This is achieved for a large number of low-lying levels, including those for which exact Bethe 
Ansatz (BA) data arc difficult to obtain due to strings of the BA quasimomcnta in the complex plane. Furthermore, 
we express the lattice eigenstates in terms of their bosonic counterparts. Since these states are then eigenstates of a 
free bosonic theory, our results have the potential to calculate expectation values of local operators in excited states. 

We consider the Hamiltonian of the XXZ model 

N 

H = £(S?SJ +1 + S P"+i + AS j S Ui)> (!) 

with periodic boundary conditions and N lattice sites. We restrict ourselves here to the critical regime — 1 < A < 1; 
data for the isotropic point A = 1 are given in the appendix. 

By a Jordan- Wigner transformation, Eq. (1) can be mapped to a model for itinerant spinlcss fcrmions. The 
corresponding Hamiltonian reads 
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where rij = CjCj and M is the total number of particles, i.e. the eigenvalue of the total number operator X)j=i n j> 
which commutes with Hf. Thus for an odd (even) number of particles, the boundary conditions of Eq. (2) are cyclic 
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(anticyclic) . The models in Eqs. (1), (2) have significant experimental relevance, either in crystals with a strongly 
anisotropic spin exchange 3,4 or in quasi one-dimensional itinerant electron models like carbon nanotubes. 5-7 Most 
recently, central quantities like the dynamical structure factor 8 and the local density of states 9 have been calculated 
for the lattice model (1) from sums over contributions of individual states. 

Historically, the model (1) has been studied extensively as a prototypical interacting many body quantum system. 
The exact solution for A = — 1 was found by Bethe; 10 Hulthen described the isotropic antiferromagnet A = l. 11 
This solution was generalized to arbitrary A by des Cloizeaux and Gaudin. 12 From these works, the ground state and 
the ground state energy were derived by Yang and Yang. 13 Excitations above the ground state were constructed by 
Takahashi (for a review, see Rcf. [14]). 

Whereas those works rely on the exact solution of the lattice model, a field-theoretical approach revealed that 
excitations with an energy AE <C 1/N above the ground state can be described asymptotically, that is in the 
limit of large chain lengths, in terms of free quasiparticles that obey bosonic statistics (for a review, see Ref. [15]). 
The corresponding effective Hamiltonian is the Gaussian model, which leads to degeneracies between certain bosonic 
excitations. Interactions between quasiparticles are captured in additional irrelevant operators, 1 which yield nonlinear 
contributions to the spectrum and generally lift the degeneracies. 

Conformal invariance relates the finite-size scaling behavior of each eigenenergy to the scaling dimensions of the 
operators in the effective field-theoretical model. 16-18 This connection has been used to predict the scaling dimensions 
of the leading irrelevant operators. 19,20 Here, we employ finite-size scaling to demonstrate the lifting of degeneracies 
for individual levels by calculating the exact contribution of the irrelevant operators to lowest order. 

Therefore, we first concentrate on those low-lying excitation levels that can be computed numerically from the BA 
for arbitrary N without convergence problems. These are mostly parameterized by real BA quasimomcnta. After 
having done this check, we can use our method to calculate low excitation energies of the lattice model with an 
accuracy of 0(N~ 4 ) for A < 1/2, and 0(N~ 2K ), K = tt/(tt - arccosA) for 1/2 < A < 1, without using the BA 
equations, regardless of the underlying BA quasimomcnta distribution. Thus especially for A < 1/2, we obtain very 
accurate analytic results for low-lying excitation energies even of relatively short chains without having to deal with 
the BA equations at all, thus also avoiding string solutions. 

Using this procedure, we then tackle the so far unanswered question of how lattice eigenstates are expressed in 
terms of bosonic states. This will help to construct the "physical" excitations seen in ab initio numerical methods or 
experiments as linear combinations of bosons. 9 

The remainder of this paper is organized as follows: Section II starts with a pedagogical introduction into the A = 
model and then treats the BA solution for general A. In the third section, the effective low-energy solution from 
bosonization is presented, including leading and higher order contributions. The lattice eigenstates and eigencnergies 
are expressed asymptotically through the eigenstates and eigenenergies of the bosonic low-energy effective Hamiltonian. 
The relative error in this asymptotic expansion is illustrated in section IV from a numerical finite size-analysis. 
Appendix A illustrates the conformal towers at the special points A = 0, 1/2, 1, and appendix B contains a table that 
illustrates our labelling of the BA levels for N = 8. 



This section summarizes the exact solution of the lattice model. In order to introduce our labeling of the energy 
levels in the large- TV-limit, we start with a pedagogical introduction into the non-interacting case A = 0. The BA 
solution is presented afterwards. 



II. THE EXACT LATTICE SOLUTION 



A. 



The free model 



After a Fourier transformation of Eq. (2) , the energy levels for the free model are given by 




(3) 
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where the M-many phases rij are chosen out of the N possible values — (N + l)/2 + p, p = 1, . . . , TV. We denote 
the values which are not occupied by a phase as "vacancies". Here, we restrict ourselves to even TV, where S z = 
0, 1, . . . , TV/2. The lowest state in each S^-sector (or equivalently, with M many particles) is given by the dense and 



FIG. 1: Sketch of the ground state a) and of a spin excitation S z = 1 in b), where one phase is removed and the remaining ones 
are again grouped symmetrically without vacancies. Here, the ground state a) is composed of an even number of particles, the 
excited state b) consists of an odd particle number. 




L R 



FIG. 2: Sketch of a current excitation m = 1. The shift of all phases can be either to the right, as depicted, or to the left. The 
left (right) Fermi points are denoted by L (R), respectively. 



symmetric distribution of the phases, without any vacancies, see Fig. la). This means that the rij in Eq. (4) are 
integer (half- integer) for M odd (even), corresponding to cyclic (anticyclic) boundary conditions in Eq. (2). The 
ground state energy Eq has no net magnetization, i.e. M = N/2, and its leading terms in a large- iV-expansion are 
then given by 
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0(N~ 5 ) . (5) 



Since we want to focus on the low-lying excitations above the ground state, it is convenient to introduce a different 
labeling of the energies. Instead of labeling each energy by the whole set of momenta as in Eq. (3), we introduce 
a notation to distinguish between three different types of excitations, according to the distribution of the rij with 
respect to the ground state. These labels refer to the true lattice states; their relation to the quantum numbers of the 
bosonic effective theory will be given in section III. 



Spin excitations: The M outermost phases are removed and the remaining ones are again placed symmetrically 
and without vacancies around the origin, Fig. 1. We will distinguish these excitations by the total spin S z = 



N 
2 

phases 



M. Due to spin-flip symmetry, the addition of M phases is energetically degenerate to the removal of M 



• Current excitations: The whole set of numbers is shifted to the right or left by m integers, Fig. 2. We label 
these excitations by m, which we take to be positive (negative) if the shift is to the right (left). 

• Particle-hole excitations on top of the "zero mode" spin and current excitations: These can always be described 
by a "shift" of occupied states relative to the filled Fermi sea, 21 labelled by two sets of integers {m!^}, {m^}, 
where m^' R denotes the number of Fermions that are shifted by n phases at the left (L) or right (R) Fermi 
point. A particular example is given in Fig. 3. These integers resemble bosonic occupation numbers but are 
so far used as unique labels for lattice eigenstates. The connection to the labels in the continuum bosonic field 
theory will be established in section III. 

In the following, we will label the energy levels by the set of numbers (S z , m, {m^} U {to^}). 

In the presence of "zero mode" excitations only, the energy of the state with 5 z -many particles removed from the 
ground state and of shifting all the phases by an integer m is obtained as 
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FIG. 3: Sketch of a particle-hole excitation labeled by (if, 2f ). The arrows denote possible shifts in the position of the phases 
in order to create excitations. The excitations do not have to be symmetric: Allowed are either shifts at the right or the left 
Fermi point or at both of them. The left (right) Fermi points are denoted by L (R), respectively. 



We now look at the particle-hole excitations, first with S z = m = 0, and we assume that at the left Fermi point, 
Nl many phases are involved in the excitations. Let j = denote the leftmost phase in the ground state. Then the 
pattern of excitations is described by the set {ji, . . . ,jN L }, where each integer takes a value ji < Nl — 1, including 
negative integers, and no two integers are equal. The analogous construction holds for the right Fermi point. The 
corresponding energy reads 

£o+ E^n^Tr- E S mV±±n + (L~R). (7) 

{r ,xj 

By considering explicit excitations on a linear spectrum, it can be shown that in leading order in I/TV one obtains 21 

27T °° 

E(0, 0, {m£} U {m«}) - E = — £ n(m L n + m*) + 0(N~ 3 ) (8) 

n=l 

Corrections of 0(N~ 3 ) cannot be written in terms of {m^} U {m^} only, as will become clear in section III. Most 
importantly, one sees that the above contributions (6), (8) can be combined linearly in leading order, 



E(S*,m, {m L n } U {m«}) - £ = _ ^ + ™ 2 + V n{m L n + m«) + 0(N- 3 ), (9) 




and the higher order terms 0(N~ 3 ) contain non-linear contributions where the excitations mix. These terms will be 
included on a more general footing in section III. 

The particle-hole contribution in Eq. (9) can be considered as the eigenvalue to the Hamiltonian 
J2 n n ( a n^ a n + a n tfl n) with bosonic operators [a^,a^J] = 6 n>m S VlfJ ,, where v, /i can be R 7 L. 15 ' 21 However, due to 
the higher-order corrections in (5), (6), (9), the labels (S z ,m,{m^} U {m^}) introduced in this section arc not the 
conventional bosonic occupation numbers. Hence the open question arises what the linear combination of the bosonic 
eigenstates is that yields the original particle-hole eigenstates of the lattice model. 

In section III we will show that such mixings of bosonic states can be determined by taking into account higher- 
order corrections to (9) and of finite interactions, A ^ 0, if the one or other lead to a splitting of the corresponding 
energies. 



B. The Bethe Ansatz solution 



In this section, the exact BA solution for the spectrum of the ArXZ-Hamiltonian with A ^ is presented. With 
its help, we demonstrate how to use the notation for the excitations introduced in the previous section. For original 
references, we refer the reader to the book by Takahashi [14]. 

The energy eigenvalues are parameterized by quasimomenta kj as 

AN M 
E ({kj}) = — -AM~J2 008 k > 

This equation looks very similar to (3). Now, however, the M = % — S z quasimomenta kj are solutions to coupled 
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algebraic equations, the Bcthe equations 



M 



(_l)M-i TT cxp [i(fc 3 - + ki)} + 1 + 2A exp [ik 3 ] 
J-l cxp [i(fcj + h)] + 1 + 2A exp [ifc ; ] ' 



(10) 



The quasimomcnta fcj pertain to interacting magnons above the ferromagnetic state. For numerical calculations it is 
more convenient to deal with the logarithmic version of these equations, 



M 



kjN = 2im.j + 2 arctan 



Asin[(fcj - k{)/2] 



cos \(kj + h)/2] + A cos [(kj - fc,)/2] 



(11) 



The lowest energy in the sector with S z = N/2 — M is given by a symmetric choice of the BA numbers in Eq. (11), 
such that 



M+l 



+ 3, 



j = l,...,M 
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In the thermodynamic limit, the ground state energy per lattice site, £o : is given by 

sinh(-7r/7 — 1) 
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We now want to obtain the energies of the lowest excitations directly from the BA equations. Some of these have 
been discussed in previous works. 2 ' 19,20 ' 22 Our aim here is to perform a systematic study starting with very small 
interactions, |A| <§C 1, and then to generalize these results to arbitrary values of A. The distribution of phases {rij} 
in Eq. (11) then defines one state uniquely. 

We want to show how to use the labels for the eigenencrgies, introduced in the previous section for A = 0, also in 
the interacting regime. As a motivation, let us first expand Eq. (11) to first order in A, i.e. close to the non-interacting 

point. This case can still be treated analytically. The corresponding quasimomenta are denoted by fcj ^ 1 ^ and are 
given by 



k 
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This expansion relies on Ac 
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<C 1 and therefore has to be taken with care for values 



±7r. One special case when this happens is near the Fermi points k 



(0) 



,(o) 



±7r/2. That the expansion 



of the BA equations in the interaction parameter cannot be trusted near the Fermi points is well known from other 
BA solvable models. 23 However, global quantities which are obtained from summing over all Bethe numbers, like the 
energy eigenvalues, turn out to be correct. 8 ' 23 More generally, the condition k^ + k^ = n defines critical pairs, 24 
with roots that can be either real or complex. The lowest excited states where these occur are current excitations 
in the S z = sector and particle-hole excitations with m^' R = l^'^, also in the S z = sector. A careful analysis 
shows 24 that these critical pairs can lead to BA numbers different from the phases that one would calculate directly 
at A = 0. We illustrate this point in appendix B for a few low-lying states in the chain with N = 8, A = 1. However, 
for any finite A, one can still label uniquely each state by the distribution of phases given by (S z , m, {m^} U {m^}) 
that one would obtain directly at A = 0, irrespective of the presence of critical pairs. This leads to general expressions 
for the energy levels to linear order in A. 

Namely, the leading terms of the spin excitation energies to linear order in A are given by 



E(S Z , 0,0) - Eo = ^(l + -A) {S z f , S* = ^ 



M. 



(15) 
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For the lowest current excitations for M = N/2 one obtains the corresponding excitation energy in linear order 

£7(0,±l,0)-£7 = ^, (16) 
which turns out to be unaffected by A in this order. Finally, the lowest particle-hole excitations have an energy 

£7(0, 0, {m L n } U {m*}) - £7 = ^ (l - V) £ n{m L n + m«). (17) 

From Eqs. (15)-(17) it is clear that A^O generally lifts the degeneracy between the lowest particle-hole and current 
levels. 

III. THE BOSONIZATION SOLUTION 

In this section, we first review the leading order of the effective bosonic Hamiltonian for the low-energy excitations 
which is accurate within 0(N~ 2 ). In the second part, next-leading corrections are included. The aim of this section 
is as follows: The eigenstates of the lattice model, \S z ,m, {m^} U are labelled by phase configurations, as 

described above. On the other hand, as will be made clear below, the eigenstates of the effective model, \S Z , m, {m^ } U 
{to^})b, arc labelled by the "zero modes" and bosonic occupation numbers as derived in [15,25] and shown in Eq. (22) 
below. Here, we wish to find those linear combinations of the bosonic states that yield the lattice eigenstates. 

A. The leading order: Non-interacting excitations 

Using conventional bosonization, the leading contribution to an effective Hamiltonian, together with its eigenenergies 
and eigenstates, for the low-energy excitations of Eq. (1) has been derived. 15,25 This Gaussian model reads 

N 1 

- ^(l + l) + £"( a « ta « +fl « ta «) ™ 

with eigenenergies 

AE (S Z , m, {m L n } U {m«}) = [E(S z ,m, {m L n } U {m*}) - Ne n ] + ^ (20) 

= -((S 2 )7A' + A-™ 2 )+5>K+m«), (21) 

Z n=l 

where So is the ground state energy per lattice site, given in Eq. (13). The effective Hamiltonian and the energies carry 
an index o to indicate that they are the leading order in an asymptotic expansion for large N and small AE <C 1/N. 
The eigenstates to Eq. (19) are given by 



|^,m,{^}U{<}) B = e iVW^« + V^ mVB ) n {aLJ yn n (afr „ |Q); (22) 



n=l 

with the following commutation relations 



[(p ,<po] = -U Q,$a =i; n, <^o =i; [cC<|] = <5n,m<W) (23) 

where fi, v stand for the superscripts R, L. 

The exponential in (22) creates the "zero mode" excitations, labelled by the integers S z ,m. The product over 
bosonic operators in (22) creates bosonic excitations, where the numbers m^' R are the bosonic occupation numbers 
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of the n-th level. The constants used in Eqs. (21), (22) are 



v 

K 



7T VT 



2 arccosA 

7T 



arccos A ' 



(24) 
(25) 



For weak interactions, v = 1 + 2A/-7T + 0(A 2 ) and K = 2 — 4A/7T + C(A 2 ), which, together with Eq. (21) agrees with 
Eqs. (15)-(17). 

In this context, one should note again that in Eqs. (20), (21) we used the same symbols as in Eqs. (6), (8) and 
Eqs. (15)-(17) by which - in the asymptotical regime - we already identified those energies from the exact solution 
with the ones from bosonization. However, as stated above, the symbols have different meanings: For the lattice 
eigenstates, they encode the phase configurations, whereas for the bosonic states, they encode bosonic occupation 
numbers. In the asymptotical regime, the exact BA eigenstates are linear combinations of the states (22) in the 
degenerate subspaces. This will be made explicit in the following section. 



B. Lifting of degeneracies due to irrelevant operators 



The Hamiltonian Eq. (19) constitutes the leading order in the large iV-limit. In the bosonization procedure, it 
results from taking account of spin-density- and spin-current-fluctuations above the ground state, where forward and 
backward scattering are included. However, Umklapp scattering processes have been neglected so far. Furthermore, 
Eq. (19) relies on the linear dispersion approximation of excitations. 

Umklapp scattering and non-linear effects in the dispersion relation induce additional terms in the low-energy 
effective Hamiltonian. These terms are expressed through bosonic fields 



4>r( x ) = 0R.O + Q R - + ^2 

1 n=l 

oo 

<j> L (x) = <Pl,o + Qlj + ^2 
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(26) 

(27) 
(28) 

In the following, leading and nextleading Umklapp processes are encoded in operators H{ ; , leading band curvature 

(29) 



(p{x) = 4> R (x) + 4>l{x). 

The operators encountered in Eqs. (19,23) are given by Q = Qr+Ql, n = Qr—Ql, = <f'Rfl+4>L,o, = 4>r. 



effects are captured by an operator H r }' 2& In our notation, these operators read 



N J \ v cos 8ttKi/(p(x)^ dx 



12N 2 ) 



+X- 



3- (1/K + K) 
2ir 



{dl4> R f : +(L ~ R) 



1 



(30) 



where v is a positive integer. Operators in Eq. (30) with :: are normally ordered. For the leading operators H r , 
the constants Ai, A+, A_ are known, 1 



Ai 

A, 



T(K) ( T(l + 1/(2^-2)) 
2tt 2 \2y/^T(l + K/(2K- 2)) / 

1 r(3 J ftT/(2 J ft:-2))r 3 (l/(2A"-2)) 
12ttA T(3/(2A - 2)) r 3 (A7(2A" - 2)) 

— tan [7rK/(2K - 2)] . 

2-7T 



(31) 
(32) 
(33) 
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FIG. 4: (Color online) The exponents 2,4 (horizontal black lines), 2K — 2 (red), AK — 4 (green dashed), 6K - 6 and 8K — 8 
(blue dotted and dotted-dashed, respectively). The vertical dashed black lines are at the values A = cos J, cos ^, cos -|, cos j. 
The open circles denote the points where logarithmic contributions occur. 



The constant Ai is given with respect to the CFT-normalization, 



lim l ^ v ( x ) c -i av (y)y = 1 (34 ) 

N-too \x — y\ za 



where 



a 2 



d = - (35) 

is the scaling dimension of the operator c lav . The expectation value in Eq. (34) is taken in the ground state. 

In the following, we will calculate the contribution of the operators (29), (30) in first order perturbation theory. 
Before going into the details of the calculation, let us first discuss the different contributions that are to be expected 
from a perturbational treatment. 



1. Scaling dimensions and perturbation theory 

The scaling dimension d that governs the behavior of the operators (19), (29), (30) under RG-transformations can 
be read off from the exponent of two-point correlation functions like (34). The fixed-point Hamiltonian (19) has 
scaling dimension d = 2, whereas the scaling dimension of the leading Umklapp operator (29) is d = 2K, see Eq. (35). 
The curvature-like term (30) has scaling dimension d = 4. 

In first order perturbation theory the correction to Eq. (21) from additional operators generally scales like N~( d ~ 2 ^ 
for finite N. Thus the operators (29), (30) can induce additional terms scaling like N~ 2 . N~^ 2K ~ 2 \ i respectively. 

On top of these leading contributions, next-leading terms exist. On the one hand, these stem from second order 
perturbation theory, giving rise, for example, to terms scaling like N~ 4 , N~ < -- 4K ~ 4: \ On the other hand, higher order 
operators can contribute in first order perturbation theory. For example, the second order Umklapp operator can yield 
a term ~ ]\[-( SK - 2 ). Generally, any additional operator can create terms in any perturbational order. We illustrate 
the exponents of the first few leading contributions in Fig. 4. 

Note that at special values of A, the exponents cross. At the free fermion point A = the amplitudes of the 
Umklapp operators vanish. However, a non-trivial crossover happens at A = 1/2. Here, the leading exponent of 
second order perturbation theory in H c , AK — 4, crosses with the exponent 2 stemming from H r . At this crossover 
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point, the two algebraic corrections merge to form a logarithmic contribution. 1:27 The same happens for higher orders. 
At the crossover points that are given by roots of unity, A = cos^-, n > 2 integer, non-trivial degeneracies between 
excited levels persist. At these points, the XXZ-mode\ has an invariance under the loop algebra sh, leading to 
additional degeneracies. 28 An unambiguous treatment of these special points has been derived within the BA. 29-31 
In appendix A, we illustrate how these degeneracies show up for the eigenvalues of the Gaussian model, Eq. (21), by 
sketching the conformal towers at A = 0, 1/2, 1. 

We shortly comment on the isotropic case. Obviously, the distinction between leading and next-leading corrections 
from Umklapp operators does not make sense for the isotropic point, A = 1. Here, all exponents 2vK — 2u — > 0, and 
the scaling dimension of He is 2 at that point. RG-studies 1 ' 2 ' 32 have shown that the corrections AE\ to the levels 
AE read 

A^ (0,0,0) - AE (0,0,0) = ^— , (36) 

V ' K ' 81n 3 A ' 

AE 1 (0 ! ±1,0)-AS (0,±1,0) = (37) 

In N 

A£i (±1,0,0)- A£ (±1,0,0) = (38) 

The amplitude g\ is known exactly, for the amplitudes g±, numerical calculations were performed. 2 

In the following, we will concentrate on the first-order contributions of the operators (29), (30). Especially, we will 

show that the operator H^ leads to the symmetric/antisymmetric combination of states with ±|to| and lifts their 
degeneracy in first order perturbation theory if S z =0. The corrections to Eq. (21) are then of order N 2 ~ 2Km ~ . In 
all other cases, this operator only contributes in second order perturbation theory, yielding corrections to Eq. (21) of 
order N A - AKm . The operator (30) always contributes in first order perturbation theory, resulting into corrections 
~N- 2 . 



2. First-order contributions from Umklapp- operators 

Let us first consider the "zero mode" states \S Z ,m,m^ = = 0)b,l, constructed according to Eq. (22). By 

inserting the mode expansions (26)-(28) into Eq. (29) for the leading Umklapp operator He , one calculates the 
expectation values of this operator between the "zero-mode" excited states. Using the commutation relations (23) 
one can show that 

B (0,l,0|ffW|0,l,0) B = 0, B <0,l,0|ff( 1 )|0,-1,0) B = C, (39) 
b(0,-1,0|^ 1) |0,-1,0) b = 0, B (0,-l,0| J ff( 1) |0,l,0) B =C, (40) 

with 

C= -(2tt) 2K A 1 A 1 - 2K . (41) 
The corresponding eigenstates in this order are 

|0,1,0) L = -L(|0,1,0) b + |0,-1,0)b) 

|0,-1,0) L = -^(|0,1,0) b -|0,-1,0)b) , 

where the labels of the lattice eigenstates have been chosen according to the discussions in section II B and [24] . Thus 

AE 1 (0,1,0) = -(2 7 r) 2K A 1 7V 1 - 2A ' (42) 
A#i (0,-1,0) = (2tt) 2K XiN 1 - 2 * . (43) 

Consequently, the leading operator describing Umklapp scattering lifts the degeneracy between the \m\ = 1 states 
for S z = and = = 0, such that the symmetric combination is energetically lower than the antisymmetric 
combination. 

Within the BA, the interacting quasiparticles above the antiferromagnetic ground state are spinons. 22 ' 33 Comparing 
the symmetric and antisymmetric current excitations with the energies of the lowest S z = two-spinon states from 
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the BA, 2,24 wc conclude that the symmetric (antisymmetric) combination of current excitations corresponds to the 
lowest two-spinon triplet (singlet) state with S z = 0. 

Particle-hole excitations at \m\ = 1 can be included as well. To determine the expectation values of the Umklapp 

operator He in Eq. (29) between these states, one again uses Eq. (22) together with the commutation relations (23). 
This results in 

|0, 1, if ' L )l = ^= (|0, 1, if L ) B + |0, -1, if - L ) B ) (44) 

|0, -1, if ' L )l = ^= ((|0, 1, if ' L ) B - |0, -1, if - L ) B ) , (45) 

with energy contributions 

A£i(0,±l,lf' L ) = (l-2A')A£i(0,±l,0) (46) 

A^(0,±l,lf if) = (1-2A-) 2 A£ 1 (0,±1,0). (47) 

The lattice eigenstates are again the symmetric and antisymmetric combinations of the bosonic eigenstates in this 
order. 

Let us now consider the states |1, ±1, 0)b, that is, states with one spin- and one current-like excitation. Proceeding 
similarly as we did in order to arrive at Eqs. (42), (43), but now including the additional spin excitation, we find 

B (l, 1,0|#«|1, 1,0) B = 0, b(1,1,0|^ 1) |1,-1,0) b = 0, (48) 
B(l 1 -l 1 0| J ff«|l ) -l ) 0) B = 0, B (l,-l,0|ff«|l,l,0) B = 0. (49) 

Especially, the terms in the second equations in (48), (49) now vanish due to the finite magnetization. This argument 
can be generalized to arbitrary S z ^ 0. Thus we conclude that for states carrying current-like |m| = 1 and spin- like 
excitations, the Umklapp operators contribute in second order perturbation theory only. 

The same is true for states without "zero mode", but bosonic excitations only. If these states are degenerate 
with respect to the Hamiltonian (19), these degeneracies are not lifted by in first order perturbation theory. 
However, second order perturbation theory generally yields a finite contribution and can thus lead to a lifting of those 
degeneracies. 

3. First-order contributions from curvature-like terms 

The operator H r in Eq. (30) yields a finite contribution for all states in first order perturbation theory. It will 
generally split the degenerate bosonic levels with the same excitation energy. In particular, we obtain for the lowest 
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levels by a straightforward evaluation of the expectation values 
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7T 
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A^i(0,0,2f' L ) 
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-6(1 
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720 5 
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(i 



AE 1 (1,0,1^) 



A£i(l,0,l£ ,L ) 



AE 1 (l,0 ) 2f'- L ) 

A^a.o.ifif) 



^2 

ASi (1,0,0) 

_2 



(60 - 60AT + (5 + 46)A" 2 )A_ 



- - (1 - 46) ) A_ 



AT 



A J B 1 (1,0,0) + — [(1-6/JQA+ 

+ (l8y/8K+{l - 2b) 2 K 2 - (36 + 12K + 60bK)) A_ 



AS 1 (1,0,0) + — [(1-6/JQA+ 

I8y/8K+(1 - 2b) 2 K 2 + (36 + 12K + 60WQ) A 
AE 1 (1,0,0) 



AT 



A r ~ 



3A' 



[(6 + 11A")A+ + (36 + 6(46 - 1)AT)A_] N~ 



(50) 
(51) 
(52) 
(53) 
(54) 

(55) 
(56) 

(57) 

(58) 
(59) 



Here we defined b := 3 — (1/AT + K). In this order of A" 2 , with 1/2 > A ^ fixed, the corresponding eigenstates are 

(60) 
(61) 
(62) 
(63) 
(64) 
(65) 
(66) 
(67) 

Thus for the above low-lying levels, the lattice eigenstates are just the bosonic eigenstates. However, the bosonic 
states have to be combined appropriately to yield the correct lattice eigenstate for the following levels: 



|0,0,0) L = 


|0,0,0) B 




|0,0,lf' L ) L = 


|0,0,lf' L ) 


B 


|0,0,lf L ) L = 


|0,0,1^ L ) 


B 


|0,0,2f' L ) L = 


|0,0,2f' L ) 


B 


|0,0,lflf) L = 


|0,0,lfl[ 


>B 


|1,0,0) L = 


|1,0,0) B 




|l,0,lf' L ) L = 


IM,lf L ) 


B 


|l,0,lflf>L = 


|1,0, lflf 


>B 



|l,0,lf L ) L = cosa|0,l,lf L ) B +sina|0,l,2f L ) B 
|l,0,2f' L ) L = -sina|0,l,l*' L ) B +cosa|0,l,2f' L ) B 



tana 



y/(26-l) 2 + 8/7i +(26-1) 
y/(2b - l) 2 + 8/K- (26-1)' 



(68) 
(69) 

(70) 
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FIG. 5: (Color online) Mixing of bosonic states according to (68), (69). The circles denote the lattice eigenstates for A = 
0, 0.1, . . . , 0.4. The red dots are the symmetric and antisymmetric combinations for A = 0. 



In complete analogy, the effect of H r in Eq. (30) on the current-carrying states with \m\ — 1 be treated. For example, 



£a(l,±l,0) = ± 



2- |) ( 3 -|) 



-Jr(60 - 60K + K 2 (365 + 4b + 60K(K - 1)))A_ 
K z 



N~ 2 (71) 
|1,±1,0) L = |1,±1,0) B . (72) 



An important consistency check of (50)-(59) is the limit A = 0. In this case, K = 2, b = 1/2, A+ — and A_ = 1/6. 
Then the above corrections yield those obtained in Eqs. (5)-(7). (Note, that the above energy correction has to be 
multiplied by 2irv/N to obtain the contribution to the total energy). The mixing of states in equations (68), (69) is 
illustrated in Fig. 5. 

In the BA solution, the energy (52) is encoded by a complex string, which makes a finite-size analysis, especially 
for large TV, difficult. The result (52) gives the leading non-universal contribution to this energy analytically, avoiding 
any problems with strings. 

Let us now look at the lattice states |0, 0, 1^'" L )l, |0,0, 2f' L ) L on the one hand and |1,0, 1^ L )l, |l,0,2f' L ) L on the 
other hand. The corresponding lattice and bosonic eigenstates are given in (62), (63) and (68), (69), respectively. For 
S* = 1, the two bosonic states |1,0, 1^' I ')b, |l,0,2f' L )B are mixed due to H r with the rotation angle given in (70). 
This angle tends to 7r/4 in the limit A — > 0, such that (68), (69) are the antisymmetric and symmetric combinations, 
respectively. A similar mixing was found from numerics for the model (1) with hard wall boundary conditions. 9 

Surprisingly, in this order, the analogous bosonic states for S z = (62), (63) do not mix for a given interaction 
1/2 > A ^ 0. The bosonization procedure directly at A = constructs fermionic states on the lattice that are the 
symmetric and antisymmetric combinations of the corresponding bosonic states. 21 However, at A = 0, |0,0, I^'^l 
and |0, 0, 2^' L )l are exactly degenerate, as can be seen from Eq. (7), such that any combination of them is allowed. 
A finite interaction lifts this degeneracy, but does not mix their bosonic equivalents in first order of H r , Hi. 

We summarize the above results for the excitation energies and give their scaling behaviors, including the leading 
pcrturbational corrections. 

ABi(0, ±1, {m£} U {m r f }) = c«(0, ±1, m£, m*)N 2 - 2K + c r (0, ±1, m£, m^N' 2 

+( f)(0,±l,m£,m«)iV 4 - 2 * (73) 
AE 1 (0,0,{m^}U{m^}) = c r (0, 0, m£, m*)N~ 2 + c^(0, 0, m£, m^)N A ~ iK (74) 
A£i(±l, m, {to^} U {m%}) = c r (±l, m, m^m^)N- 2 + (±1, to, to£, m*)N 4 - 4K , 

(75) 

where in the last equation to = 0, ±1. The terms c r , Cc are due to first order contributions from H r , Hi and have 
been calculated in Eqs. (42), (43), (46), (47), (50)-(59), (71) for the lowest states. Along that way, they can be 

(2) 

determined for any state. The terms c y c stem from second-order contributions in H c and arc not considered in this 
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FIG. 6: (Color online) Comparison of the bosonic eigenlevels (50)-(59), dotted black, with the BA values for the corresponding 
particle-hole excitations, straight red, as a function of N for A = 0.2. In the left panel, S z = 0, and S z = 1 in the right panel. 
The black dashed line in the left panel is the field-theoretical prediction for an energy given by a string-solution in the BA. 
Also shown (fat blue lines) are current excitations. The dotted-dashed line in the left panel is the field-theoretical result (43). 
The field-theoretical result (71) in the right panel cannot be distinguished from the numerical data. 



work. For the ground state at given S z , this contribution was determined in Ref. [1]. 

Therefore, the eigenvalues of the Heiscnbcrg chain can be calculated for A < 1/2 including the order j\f~ max i 4 ^ K ) 
without using the BA and thus avoiding strings completely. 



IV. NUMERICAL RESULTS 



For those states listed in appendix B, we performed a systematic finite-size scaling analysis up to lattice lengths of 
N = 2 ■ 10 3 in order to confirm our results (42)- (47), (50)-(59) and (71). 

In Fig. 6, we compare the field-theoretical results for a few energies with BA data as a function of the system length 
TV for A = 0.2. As illustrative examples, we take two current-carrying excitations, Eqs. (42), (43) and the bosonic 
states, (50)-(59), in the sectors with S z = 0, 1. 

For S z = 0, the different exponents of finite-size contributions to the energies of current-carrying excitations 
compared to states with bosonic excitations only are clearly discernible, see the left panel of Fig. 6. On the other 
hand, for S z — 1, the leading finite-size corrections to both current and bosonic excitations scale with the same 
exponent, as shown in the right panel of Fig. 6. 

For all states, the agreement is almost perfect for the longest chains studied. We also give the field-theoretical 
result for one particle-hole energy which belongs to a complex string in the BA. The field-theoretical prediction is 
independent of the actual position of roots and thus avoids any convergence problems of complex string solutions in 
the finite-size scaling analysis. 

In order to estimate the agreement quantitatively, we show the relative deviation in Fig. 7. The plots show relative 
deviations of the order of 10 -2 for N ~ 10, going down to roughly 10~ 3 for N ~ 100. This trend continues, and for 
N ~ 2000, the relative deviation is around 10 -4 for the energies considered here. 
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FIG. 7: Relative deviation of the bosonic eigenlevels (50-59), except for (52), from the BA values, for different ring lengths N. 
The large graph refers to S z = levels, whereas the inset shows S" — 1 levels. The dashed lines are guides to the eye. 



V. CONCLUSION 



We have calculated the coefficients of an asymptotic expansion in the inverse system length of a large number of 
low-lying excited energies. This calculation does not employ the BA directly, and thus avoids complex strings which 
are difficult to deal with numerically. Instead, the operators (29), (30) have to be diagonalizcd which involves the 
use of bosonic commutation relations only. This scheme is particularly useful for A < 1/2, where the energies are 
determined analytically within an accuracy including 0{N^ mayi ^' 2K ^). For the lowest ~ 50 eigenlevels that we have 
checked, no degeneracies remain that have not been already present in the lattice model. 

As a further outcome, the lattice eigenstates are expressed in terms of bosonic modes, again within the accuracy 
given above, for a fixed A. This representation of eigenstates is different from the BA representation of eigenstates. 

Applications of this approach can be manifold: Recently, there has been increasing interest in calculating the 
dynamic structure factor from a numerical solution of the BA equations. 8 ' 34,35 For other models, the BA solution is 
also being used to study dynamical quantities. 36,37 The calculations presented here to obtain the eigenenergies and 
eigenvectors could prove to be useful to obtain both numerical and even analytical results for form factors restricted 
to low excitation energies in the Heisenberg spin chain. The interest in these quantities is high, as underlined by 
the most recent work, 38 where the expectation value of the local magnetization between the ground state and a 
current-carrying state was computed. 



VI. ACKNOWLEDGMENT 



We thank F.H.L. Essler, S. Reyes and A. Struck for helpful discussions. Financial support from the Transregional 
Collaborative Research Centre SFB/TRR 49 of the Deutsche Forschungsgemeinschaft and the MATCOR school of 
excellence is gratefully acknowledged. M.B. also acknowledges financial support from the European science network 
INSTANS and hospitality at the Rudolf-Peierls-Centre for Theoretical Physics at the University of Oxford, where 
part of this work has been carried out. 



15 



APPENDIX A: CONFORM AL TOWERS FOR A = 0, 1/2, 1 



In this appendix, we illustrate the low-energy spectra (21) for A = 0, 1/2, 1 in an S z — A£Vdiagram (conformal 
tower). On the one hand, this shows the lifting of degeneracies at finite A compared to the A = 0-case. On the other 
hand, it demonstrates the occurrence of non-trivial symmetries at the special points A = 1/2 = cos ^ and A = 1. 

Each cross in a conformal tower is labeled by its coordinates in the (S z , AEq) plane. The corresponding quantum 
numbers (S z ,m, {m^} U {mf }) of the bosonic field theory are listed in the following. 

• For A = 0: 

(0,1) : (0, 0, If), (0,0, If), (0, ±1,0); 

(0, 2) : (0, 0, if If), (0, 0, 2f ), (0, 0, 2f ), (0, 0, if), (0, 0, if), (0, ±1, if), (0, ±1, if); 

(±1,1/4): (±1,0,0); 

(±1,5/4): (±1, 0, If), (±1,0, If), (±1, ±1,0); . 

(±1, 9/4) : (±1, 0, If If), (±1, 0, 2f ), (±1, 0, 2f ), (±1, 0, if), (±1, 0, if), [Ai) 

(±l,±l,lf),(±l,±l,lf); 

(±2,1): (±2,0,0); 

(±2, 2) : (±2, 0, If), (±2, 0, if), (±2, ±1, 0) 



For A = 1/2: 



(0,3/4): 

(0,1): 

(0,7/4): 

(0,2): 

(±1,1/3): 

(±1,13/12) 

(±1,4/3): 



(0,±1,0); 
(0,0, If), (0,0, If); 
(0,±l,lf),(0,±l,lf); 

(0, 0, If If), (0, 0, 2f ), (0, 0, 2f ), (0, 0, If), (0, 0, if); 
(±1,0,0); 
(±1,±1,0); 
(±1,0, If), (±1,0, If); 
(±1,25/12): (±l,±l,lf),(±l,±l,lf); 
(±2,4/3): (±2,0,0); 
(±2,25/12) : (±2, ±1,0) 



(A2) 
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FIG. 8: The lowest excitations for A = (A = 1/2) in the left (right) panel. The energies corresponding to the crosses in the 
(S z , A_Eo)-plane are given in (Al) and (A2), respectively. 
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FIG. 9: The lowest excitations for A = 1. The energies corresponding to the crosses in the (S z , A_Eo)-plane are given in (A3). 



• For A = 1: 



(0.1/2): (0,±1,0); 

(0,1): (0,0, If), (0,0, If); 

(0,3/2): (0,±l,lf),(0,±l,lf); 

(0, 2) : (0, 0, If If), (0, 0, 2f ), (0, 0, 2f ), (0, 0, if), (0, 0, if), (0, ±2, 0); 

(±1,1/2): (±1,0,0); 

(±1,1): (±1,±1,0); 

(±1,3/2): (±1,0, If), (±1,0, If); 

(±1,2): (±l,±l,lf),(±l,±l,lf); 

(±2,2): (±2,0,0) 



(A3) 



APPENDIX B: LOW ENERGY STATES FOR THE XAX-CHAIN WITH N = 8 LATTICE SITES 



In the following table, a few low-lying excited states above the ground state (note the ±S Z symmetry) are given 
for the XXX-chain with 8 lattice sites. The BA numbers are given as well as the quasimomenta and the lattice 
labels, according to section II B. For excitations with purely real quasimomenta, that is for phases \rij\ < n c with 
n c = (S z +N/2)/2, there are no convergence problems in the finite size scaling analysis because all quasimomenta are 
real. 24 Otherwise, string solutions occur, which have to be treated separately. The stars * symbolize the occurrence 
of such solutions, where a critical pair forms either a real or a complex string. These are separated off from the BA 
equations. Ref. [24] shows how to assign BA numbers to these critical pairs as well. 
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